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Abstract 



1 Introduction 



In this paper we present a state vector analysis of the generation of atomic spin squeezing by measurement 
of an optical phase shift. The frequency resolution is improved when a spin squeezed sample is used for 
spectroscopy in place of an uncorrelated sample. When light is transmitted through an atomic sample some 
photons will be scattered out of the incident beam, and this has a destructive effect on the squeezing. We 
present quantitative studies for three limiting cases: the case of a sample of atoms of size smaller than the 
, optical wavelength, the case of a large dilute sample and the case of a large dense sample. 

O 1 

S3 

In an atomic sample the population of a state \a) can be measured non-destructively by a phase shift measure- 
ment of an optical field which acts on a transition from the state \a). If the field is not too close to resonance, 
it does not drive transitions out of the state |o), but an initial state vector of the sample with a binomial 
distribution of states with varying populations n a will be modified, and the quantum mechanical uncertainty 
of the number n a is reduced. This effect has been demonstrated experimentally [Q, || , and further experiments 
have shown ||, that two separate atomic ensembles can be driven into an entangled state by measurements of 
the total phaseshifts on an optical field passing through both samples. 

A state preparation protocol that applies the outcome of quantum non-demolition (QND) measurements 
strongly relies on the fact that the measurement entails precisely the information that is applied in the changes 
of the state vector. When light interacts with atoms, apart from a phase shift of the incident field mode, also 
scattering out of the field mode occurs. The phase shift results from the interference between the incident field 
and the component of the scattered field in the incident mode, whereas scattering out of the incident mode is 
represented by components of the scattered wave function orthogonal to the incident wave function 

The scattered photons carry information about the state of the atoms which is not recorded, and therefore 
the state of the system will in general not be the one deduced from the phase shift measurements alone, but 
rather an incoherent mixture of the states that one would have determined if also the scattered photons had 
been detetected. 

The purpose of the present paper is to investigate the importance of photon scattering for the preparation of 
atomic states by phase shift measurements. In particular we shall derive criteria for the possibility to produce 
spin squeezed states. 

The paper is organized as follows. In Sec. 2, we introduce the concept of spin squeezing and some useful 
relations for the mean values and variances of spin operators. In Sec. 3, we present our model for phaseshift 
measurements, and we introduce a formalism that makes it possible to take photon scattering into account. 
In Sec. 4, we analyze the information given by the registration of phase shifts only. In Sec. 5, we present 
simulations and analytical estimates valid for a cloud which is smaller than the optical wavelength, and for 
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which photon scattering does not provide any information about the state of individual atoms. In Sec. 6, we 
consider the opposite case where the scattered photons can, in principle, be traced back to individual atoms in 
the cloud. In Sec. 7, we turn to the more complicated case of a large dense cloud, which turns out to present 
the most promising case for spin squeezing. Sec. 8 concludes the paper. 



2 Collective spin representation of an atomic sample 

The term spin squeezing originates in the treatment of two-level atoms as fictitious spin i particles, j = ^a, 
where a = (a x , er„, a z ) are the familiar 2x2 Pauli matrices in the basis of atomic states \a) and \b). For a gas 
of iV at atoms, the collective spin components 

i 

have mean values which characterize the polarization and atomic state populations of the gas and quantum 
mechanical uncertainties which characterize the population statistics. For precision in spectroscopy and in 
atomic clocks it is pertinent to have a large mean spin vector of the sample and to have as small a variance as 
possible in a spin component orthogonal to the mean spin. Assume that the mean spin points in the x-direction. 
Heisenberg's uncertainty relation states that 

AJ„AJ,>i|<J x )|, (2) 

For the state with all atoms in their respective j x = | eigenstate, the binomial distribution leads to uncertainties 
of AJ y = AJ Z = ■s/Nat/2, in accord with the previous inequality. It was shown by Wineland et al that 
if one can construct spin squeezed states which do not have the same uncertainty in the two spin components 
orthogonal to the mean spin, one may reduce the variance in a frequency measurement on N a t particles by the 
factor 

a _ N at (AJ z ) 2 

In m states were identified which for a given (J x ) have the smallest possible AJ Z . For large N at these are 
well represented by a Gaussian Ansatz for the amplitudes on states | J — M) with different eigenvalues of 
J z , in which case one obtains the approximate relation: 



1 W1 2(AJ Z 

(2J+l) 2/ ~"^ 8(AJ Z ) 2/ " ~'^ v 8(AJ Z ) 



(J*) = (J+ o)(! - 7T7 ,% ) ex P(-^7TTT2) ~ ^M-^TTTT^)- ( 4 ) 



Spin squeezed states may be produced in a number of different ways: by absorption of squeezed light ||, by 
controlled collisional interaction in Bose- Einstein condensates [Q or in a classical gas by coupling through 
a single motional degree of freedom or through an optical cavity field mode J|, . One advantage of the QND 
scheme, analyzed in the present paper, is the automatic matching of the capability to produce the state and the 
ability to detect spin squeezing, which is done by the same kind of measurement. To verify that the fluctuations 
in n a have been reduced, one has to to show that two subsequent measurements agree (to within the desired 
uncertainty). If one can produce a state with reduced number fluctuations by means of a QND measurement, 
one will also have the resolving power to make use of such reduction in a high-precision experiment. 



3 A physical setup for phase measurements. 

In Fig. 0, we illustrate a physical setup, where a beam of light enters a Stern-Gerlach interferometer which 
contains an atomic sample in one of its arms. By lenses, the field is focussed on the sample of transverse 
dimension X. For simplicity, we assume that the decomposition in plane waves of the incident photons is 
uniform for all angles 9 smaller than the focussing angle #o which, in turn, is so small that the incident field is 
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Figure 1: Configuration for a spin squeezing experiment. Atoms occupying a region in one arm of an inter- 
ferometer are illuminated by a component of an optical field, incident from the left in the figure. The phase 
shift of the light field due to interaction with the atoms in a specific internal state is registered by the different 
photocurrents in the two detectors. 
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homogeneous across the area of the atomic cloud. Let go = k9 j '(2^) denote the probability amplitude per 
unit surface for a photon to pass at the center of the mode {g^clxdy gives the probability that the photon passes 
in the area dxdy around the center. A second lens maps the field back onto the inital mode, and the second 
beam splitter of the interferometer recombines the optical beams for read out of the phase shift induced by the 
presence of atoms in the lower part of the interferometer. The atoms populate states \a) and \b), and the optical 
field couples off-resonantly the state \a) to an auxiliary atomic state, so that a phase shift on the light field is 
induced which is proportional to the population n a . 

We now take into account the scattering of the photons by the atoms. The scattering is the normal sponta- 
neous emission by the atom, given in the electric dipole approximation by well known angular distributions for 
photons of different polarization. Since our purpose is not to determine the angular distribution of scattered 
light, but rather to estimate its damaging effect on the atomic state preparation, we assume simply that every 
atom scatters photons isotropically with amplitude /. At low atomic saturation, the scattered photons are 
coherent with the incident field. 

We shall present a calculation in which the photons are scattered one at a time by the sample. The scattered 
part of the wave function of a single photon is entangled with the state of the atomic sample, since an atomic 
state with a definite sequence of atoms populating the state |a), |e) = |1 : e±,2 : e±, N a t ■ e/v at ) where e, = |o) 
or | b) leads to the photonic wave function 

V'scattj, = Mtygo (5) 

r 

where fg(Cl) depends on the state of the atomic sample. 
In the first Born approximation, is 

/e(0) = / Yl eiAtfl > 
i,ei — \a) 

where Afc is the difference between the scattered and the incident wave vectors. For 9 < 9o, /?(f2) — /?(0) = fn a . 

In the Born approximation, the flux of photons is not conserved. To remedy this problem, we write the 
angular part of the photon wave function far from the atomic sample (sum of the scattered wave function and 
of the incident wave function) in the following form which is equivalent to the first order Born approximation 
but which conserves the photon flux. 

, MQ)g for 6 > 6»o 

m) = < , — , |2 , (7) 

-i^V^-^scau) e l ~f n ^o\ for e<6o 



where 



<T S catt= I |/ ? (O) 50 | 2 dO (8) 



?>0o 

and fg is given by Eq. (^|) . 

If the atomic sample is in one of the states |e), this state is unchanged by the transmission of a photon 
through the interferometer. The scattering state of the photon, however, depends, on the argument e, and 
taking into account the unscattered component in the upper arm, we write this state in quantum notation as 

1^=^(1^ + 1^-)- (9) 

If the initial state of the cloud is in a superposition of different states |e) the joint state of the photon and 
the atoms becomes the entangled state 



|*> = £ C?\iP) p ,\e) = ]T C~ (|V> ref + \1>) g ) |e") 



(10) 
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It is at this stage, that the photodetection takes place. The photodetector 1 (resp. 2) of hte interferometer 
detects photons in the mode |1) = (e*> \tp) lei + e"** W inc )/V2 (resp. |2) = (e^ \ip) lei - inc ) / \/2) , with 

l^inc the ne ld transmitted by the lower path of the interferometer in the absence of atoms. The detection of a 
photon in one of the detectors sketched in Fig.|l| thus extracts the corresponding projection of the state vector 
\ty). This projection causes a non-continuous change of the atomic state vector amplitudes, 

C r - \ (e 10 + e-'VC 1 - °-«c«tt) e^^l*!") CV, detection 1, ( | 

C? -» § (e<* - e-*V(l - <7 acatt ) e i TF**-l»l') C e -, detection 2. 

We recall that n a and cr scatt depend on e. 

The probabilities, ttx and 7r 2 , to detect the photon in detector modes 1 and 2 are given by the squared norms 
of the vectors described by the new amplitudes after application of either projection according to (pT|). One 
may thus simulate the detection process by chosing one of the two prescriptions, update the amplitudes and 
renormalize the state vector. Detection of many photons is simulated by iterative updating of the state vector 
amplitudes. 

In addition to the detection events just described, we have identified the possibility for photons to be scattered 
into other directions Cl. The effect of such an event is of precisely the same character as the projections just 
described. If the photon is detected in the direction fi, the corresponding projection operator amounts to 
multiplying each amplitude of the initial atomic state with the corresponding scattering amplitude 

CV — ► /g{f2)CV, detection in direction 0. (12) 

The detector actually has a finite size and detects the photon in a mode with / centered around f2 but spread 
over JO. With 5ft <C 1/kX, the scattered wave function is constant over Sfl and the probability to detect a 
scattered photon within the solid angle Sfl in the direction is thus 

P(0)=£l/2M 2 |/(0)| 2 |C e | 2 ^ (13) 

and the corresponding change of atomic state vector amplitudes given by (|l^ ) . 

To simulate the scattering of photons we divide the surface of the scattering sphere in sections by longitudes 
and lattitudes, and we imagine detectors located in each section. The poles of the sphere are in the direction 
of the incident beam and the solid angle delimited by 9 < 9q correponding to the incident beam is of course 
not covered by such detectors since photons emitted in this solid angle go in the interferometer. For each 
incident photon, the probability to have a click in each detector of the sphere is computed. By adding the 
probabilities we compute a S catt and we determine the probabilities to detect the photon in the detectors 1 or 
2. The detector in which the photon is detected is chosen randomly in the simulation accordingly to all the 
calculated probabilities, and the state of the atoms is modified according to ( pi] ) or (|l2]). 

4 Information given by the interferometer 

We shall now analyze the states resulting from the interaction with the field and the detection of the photons. 
In this section we neglect photon scattering, and we study only the effect of photons measured in detectors 1 
and 2. We thus assume that er sca tt = 0, in which case we can rewrite the factors in (|i"T|), and obtain the state 
of the atoms after the detection of N\ photons in 1 and N p — N\ photons in detector 2, 

|^}=5^C ? cosU-7T0ginoJ BinU-7rflgin a J |e). (14) 

In this equation we have ignored a phase factor e™f na9 ° Np t k 1 which corresponds to a phase shift of the state \a) 
or a rotation around z in the spin language. To avoid such a rotation, one may apply an energy shift on state 
| a) or an alternative measurement scheme, where atoms in state \b) are also detected by optical phase shifts. 
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As a consequence of the photodetections, the populations of states with a definite number n a of atoms in 
state | a) are thus multiplied by the factors 

/ j \ 2iVi , j x 2(iV p -JVi) 

^N p {Ni,n a ) = cos f - TTgl-n a J sin f - 7rg 2 -n a J . (15) 

By differentiation with respect to n a , we find that ^v p (iVi, n a ) is peaked at values n° which obey, 



tan(0 + 7 r 5o 2 {n o Q )=±y^^l. (16) 

The values nP a correspond of course to atomic populations so that the probability for the photons to be detected 
in modes 1 and 2 after the interaction are in agreement with the ratios N\/N p and N 2 /N p observed by the 
measurement. 

To estimate the width of J 7 N p (Ni,n a ), we calculate the second derivative of !Fn p (Ni, n a ) at n Q a . We find 
that 

JL^f Np (N 1: n a )(n° a ) = -AN p Lg 2 £) F Np (N u n° a ). (17) 

So, the more photons that are transmitted, the narrower is the width of the peaks in Fn p (Ni, n a ). The width 
does not depend on the initial relative phase <f> between the two arms of the interferometer or on the result of 
the measurement. If we suppose that F is gaussian, the following equation gives the rms width of T in n a 

An 0int = 2 ) (18) 

The equation ([t6]) has several solutions due to the two possible signs, and due to the periodicity of the 
tan-funtion. If several such solutions lie within the initial binomial distribution of n a , the state obtained after 
the measurement will be a coherent superposition of spin states with different mean values of J z , i.e., a kind 
of "Schrodinger cat" . However, a change of Ni or N p — N% by unity changes the relative phase between peaks 
by 7r. Thus, with realistic photon detectors with an efficiency smaller than unity, the relative phase is unknown 
and the system is described by a statistical mixture of the states. 

In order to obtain spin squeezing, we want to ensure that only a single value of n° a inside the initial distribution 
obeys Eq.(|l6|), so that the detection unambiguosly leads to a more narrow distribution in n a . This requires 

— <§>fc 

> VN and — — > VN (19) 



Eq.(|^) shows that it is not enough to reduce the uncertainty in J z to have useful spin squeezing, one must 
also ensure that the mean value of J x remains large. The outcome of the interfcrometric detection is close to 
ideal in this respect. The resulting state vector has amplitudes on the different J z eigenstates which follow a 
Gaussian distribution very well, and the approximation (^) for the mean spin is close to the maximum possible 
value for any given variance of J z . 

Using Eqs.(4,18) we obtain: 

£ = l — T e^9of 2 N P /2k^ (2Q) 



The minimum value of £, for a fixed N at , is 



(21) 

and is obtained for the photon number N p = k 2 /^g^f 2 . This value of £ is the minimum value allowed as shown 



For a large number of atoms, the squeezing factor (|2l| ) can be really significant. The production of spin 
squeezed states by QND detection is susceptible, however, to two possible drawbacks caused by scattering of 
photons: 
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Figure 2: Results of a single simulation with 8 atoms whose spatial positions follow a gaussian law with a width 
of 10~ 2 A/27r. The angular spread of the incident beam is 8q = 0.457r ; close to the best that can be achieved, (a) 
Solid line : (J z ) as a function of the number of photons launched on the atoms. Long-dashed line : ((J z — (J z )) 2 ). 
Short-dashed : Expected evolution of {(J z — {J z )) 2 ) according to Eg. (|T^) 7 taking into account the width of the 
initial distribution, (b) Evolution of (J 2 ). The slight decrease shows that the state vector acquires only small 
components outside the symmetric subspace. (c) Evolution of the mean value of the spin in the horizontal 
direction (upper trace) and of the component of the spin along the expected direction (lower trace). 



• scattered photons carry information about J z , so that this quantity could in principle be known better 
than the width of TN p (N\,n a ), and according to Eq.(|J), the mean spin will be reduced 

• scattered photons carry information about the spatial distribution of atoms in the state \a). The state is 
then no longer symmetric under exchange of the particles, J- values smaller than N at /2 become populated, 
and the mean spin is accordingly reduced. 

We shall now turn to a quantitative analysis of these effects. 



5 Small cloud 

In this section we consider the case of a cloud of atoms confined to a region in space smaller than the optical 
wavelength. This implies that the scattered photons will not contain information about the individual atoms 
in the ensemble. They will, however, carry information about n a that is not known to the experimentalist who 
measures only the fields by the detectors 1 and 2. 

5.1 Numerical simulations 

In a numerical simulation of the detection process we place atoms randomly in space according to a gaussian 
probability distribution, and we assume an initial state where all atoms are in (|a) + \b))/y/2. No restriction is 
made on the state at later times, which is expanded on the whole space of dimension 2^"* . 

Fig.^j presents the evolution obtained for one particular history for a cloud of 8 atoms confined to a spatial 
region of dimension A/100 and interrogated by a beam which is focussed on the atoms with an angular aperture 
of 6>o = 0.457T. The variance of the distribution of atoms in \a) (i.e., ((J z — (J z )) 2 }) is plotted as well as the 
value expected from the results of section ^ (i.e., the multiplication of the initial distribution with the function 

The value of (J 2 ) keeps almost the initial maximal value of 20 which indicates that the cloud stays in a 
symmetric state. This is expected as the atoms are closer to each other than A and then it is not possible to 
discriminate between the atoms with the scattered photons. 
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Figure 3: Average over 90 histories, (a) Thin solid line: evolution of((J z — (Jz)) 2 ) as a function of the number 
of photons launched on the atoms. Note that this quantity is not measurable as (J z ) depends on the particular 
history. The dashed thin line gives the value expected from the scattering process only. It is a convolution between 
the initial distribution and a distribution of width given by Eq. (j&5[), both assumed to be gaussian. Solid fat line: 
evolution of ((J z — Jz aalc ) 2 ) where J Zcalc is deduced, for each history, from Eq. (|7^) and from the knowledge of 
the number N\ of photons detected in 1 and the number 7V 2 of photons detected in 2. This quantity agrees with 
Eq. (fJj|j shown by the dashed fat line, (b) Solid line and short dashed line : Evolution of *J (J x ) 2 + (Jy) 2 and 
(J x ). A rotation of ir/Ufg^N photons * s applied in the xy-plane to compensate for the light shift. Dashed line : 
expected behavior due to the squeezing realized by scattering. Dotted line : expected behavior if the width of the 
J z distribution were only given by the interference detection. 



The last graph plots the value of (J x ) and of sj (J x ) 2 + (Jy) 2 which is the length of the mean spin in the 
horizontal plane, and which may be larger than the x-component because of small angular spin rotations that 
occur during photon scattering. 

Figure [| shows the variance of J z and the mean value of J x and of the largest projection of the spin 
orthogonal to the z-axis. These results are obtained as the average over 90 independent realizations of our 
simulation. Both in Fig.|[, and in Fig|| we observe that the actual width in J z is smaller than the one 
concluded from the interferometric measurement : the fact that the atoms are coupled to other modes of the 
field also leads to squeezing. We recall, however, that the quantity ((J z ~ (J z )) 2 ) is not a measurable quantity 
as (J z ) depends on the particular history. To exploit the squeezing due to scattering, one has to deduce (J z ) 
for each experiment by keeping track of the scattered photons. 



5.2 Analytical estimates 

The effect of the scattered photons can be computed analytically. Taking r*j = for all i, Eq.(||) and Eq.(|l2|) 
show that the atomic state vector amplitudes are multiplied simply by the coefficient ^fa\n a after the detection 
of a scattered photon, and by the coefficient a/1 — o~in a in the absense of scattering where a\ = / 2 5o2vr(l+cos 9q) 
is the scattering probability per atom in state \a). 

After the detection of -/V sca tt out of a total number of N p photons, the wave function of the atoms becomes 

m = E °na (V^inaf— Vl - °mf*~ N " Mt |e), (22) 

n a 

so the probability distribution for n a is multiplied by 

Q(N scm , na) = (amlf^ (l - <r 1 n*)" ,, ~ JW (23) 
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This function is maximum on 



" a °— y 2tt(1 + cos 9 )g 2 pN p (24) 
and its width can be estimated by assuming a gaussian shape with the same second derivative at the peak value 

V8tt(1 + cos6o)g$f 2 N p 

This width is always smaller than the width in n a of the function T because go = k8o/(2y/ir), which is about 
the inverse of the transverse size of the beam on the cloud is always smaller than k, 

An ascat _ 1 gQ 1 1 , 2 q, 

An aiat 4 k 1 + cos(<9 ) 

In every single realization, the width of the n a distribution is thus set by the number of scattered photons. 
After averaging over the unkown number of scattered photons we recover the broader distribution determined 
by the interferometer readout Ni and N2- Since the scattered photons do not drive atoms out of or into the 
state | a), the atomic density matrix obtained by an average over the number of scattered photons has the same 
diagonal elements in the basis {| J, M)} as the pure state that one would expect without photon scattering. But, 
the coherence terms of the density matrix are different from that of the pure state, and therefore the length of 
the mean spin will be altered by the scattering events. This is seen in Fig.||b), where the three lower curves 
show the actual value of (J x ), of \J (J x ) 2 + (Jy) 2 and of the estimate (Q), based on the small variance of J z due 
to the scattering. There is excellent agreement between these curves. The upper curve shows the larger value 
of the mean spin, that one would have obtained in the absence of scattering. 

Our numerical simulations and our analytical estimates show that the effect of the scattering is to produce 
stronger squeezing than the interference measurement. After averaging over the unresolved scattering histories 
this squeezing does not affect the n a populations. But, its effect is to reduce the mean spin (J x ). In principle one 
could determine the number of scattered photons by the difference betwen the number of incident photons and 
the number of photons detected in the interferometer. For a poissonian source of light, however, this number 
cannot be determined to a higher precision than WNp, which turns out to be larger than the required precision 
on the loss in photon number due to scattering. 

As pointed out in section ||, the pertinent factor for spin squeezing is £ defined in Eq.(||). Using Eq. 
for AJ Z and Eqs.(4,25) to determine (J x ), we obtain 

£ = - e ^9lf 2 N P (27) 

VWt^Np~ng 2 £ 



The minimum value of £, for a fixed N at , is 



and is obtained for the photon number N p = l/(\/4iTg 2 f 2 ). Because the size of the incident beam is larger than 
the wave length, g < k and ^Min is larger than in the ideal case (|2l|). 



, 4e k 

?Min = \ — t- t ■ (28) 

nN at g 



6 Large dilute cloud. 

The effect of the scattering in the case of a cloud of extension larger than A is very different from the case of a 
small cloud. The number of scattered photons will still, as in the case of a small cloud, give us information on 
the total number of atoms in \a) and thus scattering will by itself produce squeezing. But in the case of a big 
cloud the angular distribution of the scattered photons gives also information on the position of the atoms in 
| a) and the atomic system therefore looses its invariance under permutation of the particles. The system will 
thus no longer be fully represented by the (N at + 1) Dicke states of maximal J = N at /2. This does not affect 
the z component of the spin but it strongly reduces the horizontal component of the spin. 
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6.1 Numerical simulation 



In order to analyse the effect of the scattering alone, we have performed simulations of the evolution of the 
atomic state under detection of scattered photons. The directions of photodetection were chosen according to 
the procedure outlined in Sec. Ill, and Fig. || presents the result of a single simulation. In this simulation, we 
see a decrease of ((J z — (J z )) 2 ) as well as a departure from the symmetric space (decrease of (J 2 )). It is also 
clear that the decrease of ((J z — (J z )) 2 } is not as rapid as suggested by Eq.(p5|) which was valid for the small 
cloud. 

Let us present a naive argument for the decrease of ((J z — (J z )) 2 ): Assume that the cloud is in a state |e) 
where each atom is either in \a) or \b). This state is unaffected by the interaction with the photon field. Let 
p denote the probability that a photon is scattered by the cloud of atoms. In the absence of superradiance we 
expect p = o\n a . The number of scattered photons N scaU thus has mean value pN p and uncertainty y/pN p , and 
we estimate n a by N scat t /crN p with an uncertainty of An a = yj n a /aN p . If the state of the cloud expands over 
different n a , the measurement of N scatt will reduce the width of the distribution by multiplying with a function 
of width An a = y/n a /aN p . Taking a = 4ng 2 f 2 and an initial state where all atoms are in (\a) + \b))/y/2 

N at + Nat 

This function is plotted in Fig.^ (a). It reproduces quite well the numerical evolution, although the numerical 
evolution shows a faster squeezing. 

After averaging over histories, as discussed in section ||, the squeezing due to scattering has no effect on the 
distribution of n a but it will reduce the horizontal component of the total spin. The upper dashed line in figure 
H(c) predicts the reduction due to this effect, according to Eq.(|) taking J = N at /2. Wc sec that the horizontal 
spin determined in our simulation is even smaller than this estimate. The reason for this is that the state of 
the cloud is not in the symmetric subspace as assumed when we put J = N at /2 in Eq.(|j). The state of the 
cloud may be expanded on subspaces with different total J, and the mean value of J x is averaged over these 
different subspace components. To check the consistency of this picture we compare the decrease of (J 2 ) with 
that of (J x )- The maximum (J) for a given (J 2 ) = (J (J + 1)) is y/1 + 4(J(J + l))/2 - 1/2 and it is obtained if 
only the subspace with J = (J) is populated. If we assume that the J z distribution (centered around zero) is 
the same in all subspaces we estimate that 

ws .-wte(-i + ^±SE±S) (30) 

The curve corresponding to the right-hand side of this inequality is plotted in Fig.||. It reproduces quite well 
the decrease of (J x )- 

6.2 Analytical estimates 

If the atomic cloud is dilute enough and not too large it is possible to know by which atom any photon has been 
emitted. In other words, it is possible to design an optical system which collects all the photons emitted outside 
the incident mode and which produces a separate image of each atom. This implies that the flux of scattered 
photons is only proportional to n a and not to n 2 : there is no superradiance. Then the scattered photons give 
information on the total number of atoms in \a) but also on which atoms arc in \a). The second effect damages 
the squeezing realised by the interference measurement as it decreases correlations between atoms. When a 
photon is transmitted, the probability that it is scattered by an atom in \a) is approximately Ang^f 2 . Thus 
after about N p = l/Airg 2 ,/ 2 photons have been transmitted, we know almost with certainty if the atom is in 
1 6) or |o), and there is then almost no correlation between the atoms and hence no squeezing. Indeed, as each 
atom is either in \a) or \b), (j x ) — for each atom. In order to observe squeezing the number of photons used 
in the experiment is limited and this in turn limits the reduction in Var(J z ). 
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Figure 4: The cloud contains 8 atoms spatially distributed according to a 3D gaussian probability law of rms width 
10A/27T. The incident beam, focussed on the atoms, has an angular spread 8q = 2^/ngo/k = 2y / 7r * 0.05 = 0.14. 
(a) : Evolution of (J z ) (solid line) and ((J z — (J z )) 2 ) (long dashed line) as functions of the number of photons 
which are launched on the atoms. The dotted line indicates the value of ((J z — (J z )) 2 } expected for a dense cloud 
of 8 atoms illuminated by the same beam (Eq. (jjqjj. The short dashed line gives the result of Eg. (jl^j. (b) : 
Evolution of (J 2 ). The first drop from the maximal value 20 occurs at the first detection of a scattered photon, 
(c) : Evolution of the mean horizontal spin (solid line). In dashed lines is shown the expected value of (J x ) 
according to Eq.^) for J = N at /2. The result of Eq. (p?Z[) is shown with the short-dashed line. 



Within a more quantitative analysis of the effect of scattering let <f> denote the flux of incident photons. As 
long as an atom scatters no photons, its state evolves according to the effective Hamiltonian 

H nh = - l -^glf 2 \a){a\. (31) 

The probability to scatter a photon during 8t is 

PSt = (a\ ip \a) 2 <j)ATtglf 2 5t. (32) 

After averaging over histories only the coherence term a ao of the density matrix evolves and it obeys 

dcr a b <t> . 2f 2 rtQ \ 
__ = --4irg Q f a ah . (33) 

Thus, (j x ) = l/2(<7 & + a* b ) follows the same exponential decay, and it follows that the total spin J x — ^j Xi 
decreases as 

(Jx)t = (J x )oe~ 2 ^ f2t = (J x )oe~ 2N ^ 2 ° f2 (34) 

The reduction of (J x ) due to the interferometric dectection, estimated using Eq.([|) is less important than the 
reduction of (J x ) due to scattering and we will neglect it here. Var(J z ) is given by Eq. (fl8|) , and the squeezing 
factor writes 

e=^^ T e 2N ^of\ (35) 

Its minimum value is 

V g 

and it is obtained for N p = 1/ '(Ang^f 2 ). 

Although the physical regimes are very different, the result is similar to Eq.(26). 
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7 Large dense cloud 



We now turn to an analysis of the situation of many atoms in a large cloud with a density which is large 
compared to 1/A 3 . In this case, we are not able to carry out simulations, and we arc therefore restricted to an 
analytical approach. We shall make the assumption of a dense cloud, so that it may be divided into a large 
number of cells of size a 3 smaller than A 3 but still containing a large number of atoms. For simplicity, we 
assume that the distribution of atoms is uniform over a box of size L x = L y , L z . As the size of the cell is smaller 
than A 3 , the scattering of photons does not bring information on the repartitioning of the atoms in \a) inside 
each cell and the state of the cloud can be described by states that are symmetric under exchange of particles 
inside each cell. The states \n) = \n\,ri2, ...,tim) with well defined number n\, i%2,...,nM of atoms in \a) in the 
cells 1,2, ...,M are not modified by the scattering. The number of atoms per cell is sufficiently large so that we 
consider that the number of atoms in \a), on the order of N ce u/2, can be considered as a continuous variable 
with fluctuations of order y/N ce u/2. 

7.1 Initial state of the cloud 

Initially, all the atoms are in the state l/^/2(\a) + \b)). The expansion of the state of the cloud on the basis 
{\n}} is then 

\tp) at = J dm J dm... j dnMc(n 1 )c(n 2 )...c(n M )\ni,n2, ...,n M ) (37) 

where c(n) is the square root of the binomial distribution of mean value N ce u/2 which we will approximate by 
a gaussian distribution. 



We will now define the new variables 



= V h Er«rcos(^.f r ) 

N s W *=\flY.Tn r sHK^r) ( 38 ) 

with 

ff= al x x° + alyip + alzZ® 

where i z < or (i z = and i x < 0) or (i z = i x = and i y < 0). 

Note that all the operators N a and Np commute since they are all diagonal in the basis 
The initial state of the cloud can be expanded on the new basis of eigenstates 

Wat= JdN fdN Cl fdN Sl ...fdN C(M _ 1)/2 fdN S{M _ 1>/2 (40) 
C(N )h(N c ^h(N s J...h(N C(M _ 1)/2 )h(N S(M _ 1)/2 )\^ 

where h is a gaussian centered on with an rms width of N ce u /2 (h 2 has a width y/N ce u/2 equal to the width 
of c 2 (ni)) and C has the same width but is centered on \[MN ce uj2. Initially, there is no correlation between 
the distributions in N s / C . and they are all the same. 

The vector space can be seen as a tensor product of (M+l)/2 spaces: The space acted upon by No, and for 
each i, the space £i acted upon by the operators N Ci and N Si . All these spaces have infinite dimension and they 
admit as basis states \N Ci ,N Si ) where N Ci ,N Si are real. This basis will turn out to be useful for the description 
of the state vector dynamics due to photon scattering. 
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Figure 5: Inital distribution among the states (1^1, ?i2> TIm)) or the states 

(\Nq, N Cl , N Sl , N C(M _ 1)/2 , N S(M _ 1)/2 ^) in the case where there are only three cells (M — 3). Repre- 
sented is a surface on which the population is constant. The circle surrounding the origin is the projection of 
the sphere on the (N c , N s ) plane. 

7.2 Effect of the scattering 

Due to energy and momentum conservation in the scattering process, a fluctuation in the cloud will couple to 
the light only if its Fourier transform has components Ki on the scattering sphere defined as the vectors kd 
which fulfill | kd — ki nc \ = k. 

We will now assume for simplicity that the Fourier transform associated with each Ki are uniform over 
disjoint volumes 5k x — 2n/L x , 5k y — 2n/L y and Sk z = 2n/L z , shown as small rectangles in Fig. |[ The circle 
in Fig. ^ depicts the scattering sphere and crosses represent the dicrete Ki wave vectors which contribute to the 
scattering. With this approximation, the observation of a scattered photon gives information on the modulation 
of the atomic population in state \a) with the corresponding discrete wave vector. 

Let us consider a Ki which is on the scattering sphere and the associated solid angle 5fli. The effect on the 
state of the system associated with the detection of a photon in this solid angle is given by the operator 

V(Q) = VSQ igo f n/^-^r = ^m igo fJK( N ^ + iNa j. (41) 
f 

This operator changes the relative phase between the component with fluctuations in cosinus and sinus com- 
ponents. Because the operator in Eq.(40) acts only on the space Si, no correlations between fluctuations at 
different wave vector appear and the state of the cloud will thus stay on the form 

|V>o) ... |V>(M-i)/ 2 ) (42) 

where \ipi) is a state of the space Ei. Note that this result is an approximation which relies on our assumption 
that the Fourier transform of the fluctuations at different Ki are disjoint. Furthermore, the assumption that 
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Figure 6: Diffusion sphere and components of the discret Fourrier decomposition which participate to the diffu- 
sion 



atoms are uniformly distributed over our spatial grid is important. Indeed, if the number of atoms in \a) 
in different cells is not the same for each cell, the initial state would present correlations between Kfs and 
knowledge about the state in the subspace would also implies knowledge about the state in other £j. 

If a photon is detected in the incident mode, the effect on the state of the atoms, to second order in /, is 
given by 

V = U\J l -9of 2 Y ^ ( N l+^)- ( 43 ) 

i 

Thus, in this approximation, neither scattering nor absence of scattering yields correlations between different 
Ki's, and the state of the cloud will stay in a product state as in Eq.(42). 

After transmission of N p photons and the detection of N sc photons scattered in the direction = 
kmc ~ Ki, the wave function in the space Ei becomes 



dN Cz J dN Si (N Ci + iN Si f*° (y 1 - g$f 2 Y SQ i ( N l + N l)J h(N Cz )h(N Sl ) \N Ct ,N Si ) (44) 

Thus, the population of the states \N Ci , N Si ) is multiplied by the factor 



G^Nl+Nl) = (Nl + Nl) N - (1 - g 2 ff SQi (J\£ + iV s 2 J) Np ~ Nsc , (45) 

which depends only on AT?. + N%. . This is expected as the detection of the scattered photons does not reveal 
any information about the phase of the spatial grating in the cloud of atoms in \a). 
A calculation similar to the one in Sec. 6.2 shows that Q has a width 

AN, = 1 (46) 

^5^fglpN p 

Figure ^ depicts the final distribution over the states Nq, (N Ci , N Si ) in the simple case where the cloud is 
divided into only three cells. 
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Figure 7: Distribution over the states (\ni,ri2, ■■■,tim)) or the states (\Nq, N Cl , N Sl , N C(M _ 1)/2 , N S(M _ 1)/2 )) 
after the measurement. Compare with Fig. 5. The case of only 3 cells is represented, the distribution of the 
population in the basis (\N C , N s ) ) for a well defined Nq does not depend on No and is given by the projection of 
the torus on the plane (N c ,N a ). 



7.3 Effect of scattering on < J x > 

Photons scattered in the forward direction in the solid angle 9 < 9q give information on Nq and thus on the total 
number of atoms in \a). The other scattered photons give information on the spatial fluctuations of the atoms 
in | a) and, in the case we consider where the number of atoms per cell is large and approximately constant, no 
information on No is given by the photons scattered outside the solid angle 8 < Bo- This is completely different 
from the case considered in section |6.2| where in a single experiment the number of atoms in \a) was mainly 
determined by the scattering. 

In the case considered here, even if No is not affected by the scattering outside the incident beam, the state 
is affected by scattering events due to the departure from the space of states which is symmetric under exchange 
of atoms. After a series of detections, the state of the cloud has components in different subspaces J because 
the scattering (or its absence) brings the system into non symmetric states with J < N/2. We thus expect 
(J x ) to be smaller than the value obtained for a state with the same distribution of J z eigenstates but in the 
symmetric subspace J = N/2. 

We have 

(■>*> =!>*«> ( 4? ) 

i 

where {J Xi ) is the mean value of J x corresponding to the atoms of the cell i. It is 

(Jx,) = J dn 1 ...dn i -idn i+ i...dn M (J dm\f(ni, ...,n M )\ 2 ^j (Jx j )ni,...,n < _i,n i+ i,...,nM ( 48 ) 
The state is invariant by exchange of atoms inside a single cell, and the results of Eq.(|j) can be used to estimate 

by 
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(T ) , ~ N cell sva^j. j c 4 q\ 

\^a;i /nl,. . . ,n^_ i ,ni+i , . . .,riK — 2 V / 

Taking i = 1, we note that Var(J Zl ) is the rms width of the function 

F( ni ) = C + n2+ ^ ?lM ) (V^M 2 + ^>i) 2 ) (50) 

where C'(Nq) is the final population distribution over the states {|ATq)} and Ti'^Ni) = Gi(Ni)h 4 (Ni) is the final 
population distribution among the states {\N C . = cos(a)N i: N s . = sin(a)iVj)} (It does not depend on the angle 

«)■ _ 

For the initial distribution, C (ni / VM +...) has an rms width of V MN C eii /2 in rti and Hi ^V^ci (™i) 2 + (^l) 2 

h 2 (y^2/M cos(Ki.fi)ni + ...)h 2 (^2/M s'm(Ki.ri)ni + ...) has an rms width of ^jMN ce u/% as a function of n\. 
Thus, Eq.([50|) gives an rms width AJ Zl equal to V N ce ii /2 for any (712, «-m), which is expected. 

After the detection, the width of C has been reduced by the interference measurement to An ainteT /vM and 
the rms width in n\ oiC'(rii/yM + ...) is 

Ani mter = An ainter . (51) 

The detection of scattered photons modifies the distribution H! i which is then no longer separable in N Ci and 
N Si . But it follows from the definition ( |38| ) and from Eq|l6] that, in average, 7^ has an rms width in n\ which 

scales as 

Ann ~ VMAiV, (52) 

If we assume that the distribution of n% is given by the initial Gaussian, which is multiplied by gaussian factors 
due to the interferometric detection and due to the scattering, we obtain the result 

~A~n{ = An 2 a . ; + ^ MANf + ^2 MN^i (53) 

i,K t eS sea tt i,Ki<t.S scatt 

The squeezing due to scattering and due to interferometric detection is significant, and we can hence ignore 
the last term, which is due to the initial width of the distribution. The contribution due to the interferometric 
detection is given by Eq.(18), and the sum over the i for which Ki is on the scattering sphere follows from 
Eq.©: 

_E m^ =2 ^ (54) 

i,Ki£S scat t 

By comparison of Eq.(18) and (|54|), we see that the contribution due to interferometric detection is approximately 
a factor g 2 , / k 2 times the one due to scattering. 1 /g 2 is larger than the area of the cloud which itself is much 
larger than A 2 . Thus, g 2 /k 2 <C 1 and the width in n\ is mainly determined by the scattered photons : 

A^f = Vak * h N r (55) 



Coming back to Eq. (|||) and Eq. (^) , we get 

2 



J x ) ~ ^e-f tit**' (56) 



Thus, the squeezing factor writes 



Its minimum value, achieved for N p = 4/ (irg 2 f 2 ), is 



£ = 7 == 2/ ef^/ iv P (57) 



6«« = a/t4^-> (58) 



which is similar to the result obtained for the small cloud and for the dilute sample. 
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8 Conclusion 



We have presented an analysis of the change of the atomic distribution on internal levels caused by a mea- 
surement of the phase shift of an optical field traversing the atomic sample. In every run of the detection 
experiment, the atomic population statistics is modified in accordance with the phase shift measurements. In 
practical spin squeezing, the reduced variance in one of the spin components is not the only relevant parameter, 
one has to observe the change of the length of the mean spin (J x ) as well. The mean spin is reduced, and in 
case of perfect detection, we estimate the optimum number of detected photons, and the optimum value of the 
squeezing parameter £ = \/l/N at , (gjj). 

In an experimental implementation, scattering of photons is inevitable. If these photons are not detected, 
they will have no average effect on the atomic populations, but scattering leads to further reduction of the mean 
spin and an increase of £. Two different mechanisms are shown to be responsible for this reduction: In case 
of the small cloud, the scattered photons carry information on n a , which forces a reduction in (J x ). In case of 
the dilute cloud, they carry information on the spatial distribution of the atoms in \a). The general problem 
is difficult to treat, and we focussed on three different limiting cases : a cloud of size smaller than the wave 
length, a dilute cloud where each scattered photon can be traced back to a single atom, and a large cloud of 
density larger than 1/A 3 . Both numerical simulations and analytical approach were carried out and we showed 
that scattering decreases the mean spin (J x ). 

Although the physics is very different, the scattering gives rise to approximately the same optimum squeezing 



which is obtained by detecting the phase shift of a given number of photons, N p ~ \/(g^f 2 ). We cannot focus 
to better than within a wavelength, hence the optimum always exceeds the ideal results ( |2l|). In the case of a 
large cloud, 1/Qq is of the order of the area of the beam at the focus, and go is limited by \/A, where A is the 
area of the cloud . In this case we can therefore rewrite the expression for ^Min 



where D — N at X 2 /A is the optical density of the cloud. Thus for a dense cloud spin squeezing is possible, 
whereas for a very dilute cloud we recover the result that the possibility to know, for any scattered photon, 
which atom it comes, effectively prevents spin squeezing. 

A better squeezing could be achieved if the atomic cloud lie in an optical cavity so that a photon passes 
effectively nt times in average through the cloud, where nt is the finesse of the cavity. In particular, this would 
enable squeezing of a dilute cloud. To estimate the best possible squeezing in this case, two opposit effects 
should be taken into account. First, the width in n a infered from the phase shift of the beam after the passage 
in the cavity is decreased by a factor nt compared to Eq.[l8| Second, the probability that a photon is scattered 
before it leaves the cavity is also enhanced by nt so that the number of photons N p that can be used before 
(J x ) decreases too much is also decreased by n t . But, as the width in n a induced by the phaseshift decreases 
only in 1 / y/N^, the best squeezing factor £ is still decreased by a factor 1 / y/rit compared to Eq.^ and could 
thus become smaller than 1. 
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